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Abstract: We compute the two-loop master integrals for non-leptonic heavy-to-heavy 
decays analytically in a recently-proposed canonical basis. Eor this genuine two-loop, 
two-scale problem we first derive a basis for the master integrals that disentangles the 
kinematics from the space-time dimension in the differential equations, and subsequently 
solve the latter in terms of iterated integrals up to weight four. The solution constitutes 
another valuable example of the finding of a canonical basis for two-loop master integrals 
that have two different internal masses, and assumes a form that is ideally suited for a sub¬ 
sequent convolution with the light-cone distribution amplitude in the framework of QCD 
factorisation. 
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1. Introduction 

Non-leptonic R-decays are interesting for a number of phenomenological applications like 
the extraction of CKM elements and the study of CP asymmetries. Their study has 
already entered the area of precision physics, both on the experimental [ 1 ] and on the the¬ 
oretical side. However, their theoretical description is complicated by the purely hadronic 
environment, entailing QCD effects from many widely separated scales. The two main 
approaches to non-leptonic R-decays are flavour symmetries of the light quarks [ 2 ] and 
factorisation frameworks such as pQCD [3] and QCD factorisation (QCDF) [4-6]. In the 
latter framework, next-to-leading order (NLO) corrections to both, heavy-to-heavy [5] and 
heavy-to-light [4,7] transitions have been known since more than a decade. More recently. 
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also next-to-next-to-leading order (NNLO) results for heavy-to-light decays have become 
available [8-12]. In the present article, we consider NNLO corrections also to the heavy-to- 
heavy decays such as B ^ Dn in the framework of QCDF [13]. In the heavy-quark limit, 
the decay amplitude for —)• D~^7t~ is given by [5] 




(I.l) 


where Ot are the operators from the effective Hamiltonian that describe the underly¬ 
ing weak decay. The form factors and the pion light-cone distribution amplitude 

(LCDA) ‘k.n-(u), with momentum fractions u and 1 — u shared among the pion constituents, 
are the non-perturbative inputs. The hard-scattering kernels Tij(u), on the other hand, 
can be evaluated in a perturbative expansion in the strong coupling, and are known in 
QCD to NLO accuracy [5]. Yet it is interesting to go beyond NLO in H —)• Dtt transitions: 
Since the contribution at NLO is colour suppressed and appears alongside small Wilson 
coefficients, the NNLO corrections may be significant in size. Moreover, since there is 
neither a colour-suppressed tree amplitude nor penguin contributions, and spectator scat¬ 
tering and weak annihilation are power-suppressed [5], we have only the vertex kernels to 
the colour-allowed tree amplitude. A precise theory prediction of this single contribution, 
together with comparison to experimental data, might give a reliable estimate of the size 
of power corrections in the QCDF framework. 

The evaluation of Feynman diagrams that contribute to the NNLO hard-scattering 
kernel amounts to the computation of ~ 70 two-loop diagrams. By using contemporary 
techniques to evaluate multi-loop integrals, the two-loop Feynman diagrams are reduced 
to a small set of a few dozens of master integrals. A powerful method to evaluate the 
latter analytically are differential equations [14-16]. This method was recently refined by 
Henn [17]. Considering that the basis of master integrals is not unique, Henn discovered 
that in a suitably chosen basis - denoted as canonical basis - the differential equations 
can be cast into a form that factorises the dependence on the kinematic variables from 
that on the number of space-time dimensions D. In this case, the solution is expressed in 
terms of iterated integrals. This method was recently applied to a number of problems for 
loop [11,18-26] and phase-space [27,28] integrals. 

To the present day, the construction of the canonical basis is mostly based on experience 
or experimentation, rather than on a systematic procedure, although developments in this 
direction have recently become available [21,29,30]. In the future it would be most desirable 
to have a general algorithm for finding a canonical basis for arbitrary external kinematics 
and numbers of loops, legs, scales, and space-time dimensions. Therefore, every non-trivial 
example of a canonical basis is most valuable, and our results contribute towards finding a 
general algorithm for constructing the canonical basis. 

Last but not least, if the master integrals that enter the hard-scattering kernels Tij{u) 
are written in terms of iterated integrals, the convolution with the pion LCDA in (1.1) 
simplifies to a large extent. Our results therefore catalyse the steps necessary to obtain the 
decay amplitudes considerably, and constitute an important step towards the phenomenol¬ 
ogy of H —5- Dtt decays at NNLO in QCDF. 
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This paper is organized as follows. In section 2 we introduce the kinematics of the 
two-body decay and present the generic form of the differential equations with respect to 
the kinematic variables. We proceed by defining Goncharov poly logarithm in section 3, 
which are a class of iterated integrals suited to describe the solutions to the differential 
equations. In section 4 the canonical basis is defined and the expressions for the master 
integrals in this basis are presented. We also elaborate on strategies to find a canonical 
basis. The boundary conditions for the integrals are discussed in section 5 and the results 
are presented in section 6. In section 7 we comment on the performed cross-checks before 
concluding in section 8. In appendix A we collect the matrices that contain all relevant 
information on the differential equations. The analytic results of all master integrals are 
also available electronically [31]. 


2. Kinematics 


We consider the kinematics of the decay —)■ which emerges from the underlying 

weak transition b —)• cud. A sample of Feynman diagrams contributing to the two-loop 
hard-scattering kernels is given in figure 1. The complete set of diagrams consists of those 
shown in figures 15 and 16 of [5], supplemented by gluon self-energy insertions in one-loop 
diagrams. All external momenta are taken to be incoming throughout this work, and 
denote the external momenta of the b and the c quark, respectively, which fulfill the on-shell 
constraints 3 = The constituents of the pion share the momentum q with gi = uq 
and g2 = (1 — 'w)'? = uq, where u G [0,1] is the momentum fraction of the quarks inside the 
pion entering eq. (1.1) in a convolution of the hard-scattering kernel with the pion LCDA. 
We consider the pion to be massless, i.e. q^ = qf 2 — the linear dependence of 

the momenta, qi + q 2 = Q = —Q 3 — Qa, the kinematics is completely determined by two of 
the on-shell conditions and one additional kinematic invariant, for instance 

2222 t/ 2 i 2 \ 

94 = > 93 = > 9394 = + '^c) ■ (2-1) 

We apply commonly used multi-loop techniques which include integration-by-parts identi¬ 
ties [32,33] and the Laporta algorithm [34], and reduce the two-loop Feynman diagrams to 
master integrals [35,36]. Furthermore, we construct the differential equation of the latter 
with respect to kinematic variables. In the derivation of eq. (1.1) the charm quark was as¬ 
sumed to be heavy. Hence, the ratio mcjmb remains fixed in the heavy-quark limit and our 
master integrals depend on two scales: the momentum fraction u and the ratio of the heavy 
quark masses 2 ; = They are further functions of the kinematic invariants (2.1) 

C(u,z) = C{u,q‘l{z),{qiq 3 ){z),ql{z),z). (2.2) 


Thus, the total derivative of a generic master integral C with respect to u is given by 


d^_(^ 

du du ’ 

whereas the one in 2 ; reads 

d^_^ c^d^ dC djqsqA) dC dgl 
dz dz (9^3 dz 9 ( 9394 ) dz dz 


(2.3) 


(2.4) 
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Figure 1: Sample of Feynman diagrams: and qs are the momenta of the quark lines with masses 

mb and TOc, respectively. qi = uq and <72 = uq are the momenta of the light quark and anti-quark, 
respectively, q = qi + q 2 is the momentum of the pion. All momenta are incoming. The black 
square denotes an operator insertion from the weak effective Hamiltonian. 


The computation of dCjdz is straightforward. The partial derivatives of C with respect 
to the kinematics on the r.h.s. of eq. (2.4) can be expressed in terms of partial derivatives 
with respect to the momenta < 73 ^^ and ^ 4 ^^ [37], which can be easily carried out. Note that 
the last term on the r.h.s. vanishes since dq\/dz = 0. We finally obtain 

dC dC If dC dC \ 

This is the differential equation with respect to valid for a generic master integral C{u, z). 


3. Iterated integrals and Goncharov polylogarithms 

The classical example of iterated integrals is given by the harmonic polylogarithms (HPLs) 
[38]. They generalise the ordinary polylogarithms and are defined by 

= f dt j (3-1) 

Jo 

where the parameters a* can be 0 or ±1, and n is the weight of the HPL. The integral (3.1) 
diverges for HPLs with trailing zeroes. In order to handle HPLs in such cases, one defines 
(x) = ln"'(x). The weight functions fai{x) are simply 

h{x) = , /o(x) = - , . (3.2) 

1 — X X 1 + X 

The HPLs fulfil a Hopf algebra according to 

Hsix)H^{x) = ^ Hsix) , (3.3) 

cGaldJfo 

where a tt) 6 are all possibilities of arranging the elements of a and b such that the internal 
order of the elements of a and b is preserved individually (cf. also [39]). Hence the product 
of two HPLs of weights wi and W 2 has weight Wi + W 2 - The Hopf algebra can also be used 
to extract singular behaviour near x = 0 or x = 1. Due to the relation 

77o,..., 0 , 1 ( 1 ) = a (3.4) 
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with k — 1 zeroes and k > 1, one also assigns the weight k to numbers like Cfe and vr^. 

A generalisation of the HPLs are the Goncharov polylogarithms [40], whose definition 
reads 


Gai,a2,...,a„ix) — 


dt 


t — ai 


Ga2,...,an (t) 


(3.5) 


and Gg^(x) = Jdj^^(x). They fulfil a Hopf algebra that has the same structure as (3.3), and 
allow for more general weights than just 0 or ±1. In particular, in multi-scale problems 
the argument x can be represented by one scale, and the remaining scales are comprised in 
the weights ai. In our problem at hand, it is most convenient to choose u as the argument 
of the Goncharov polylogarithm whenever there is a dependence on this scale, bearing in 
mind that this choice simplifies a subsequent convolution with the light-cone distribution 
amplitude, which in a Gegenbauer expansion is a n-dependent polynomial. In this case the 
weights are either integer (0, ±1) or one of the following six z-dependent weights^, 


Ol 


1 

1 - z ’ 


02 


Z-l ’ 


03 


04 


1 

1 

r+V^’ 


05 




'V^-I ’ 



(3.6) 


Goncharov polylogarithms that do not depend on u are written in terms of integer weights 
and argument 2 : or y/z. Products of Goncharov polylogarithms of the same argument are 
expanded by means of the Hopf algebra. 


4. The canonical basis 

We work in dimensional regularisation with ZD = 4 —2e and evaluate the two-loop, two-scale 
master integrals by applying the method proposed by Henn [17]. Gonsidering a specific 
power in the e-expansion of a master integral, the associated function is called uniform if 
each summand has the same weight. Moreover, a uniform function is called pure, if its 
derivative with respect to any one of its arguments yields a uniform function whose weight 
is lowered by one unit. 

The proposal in [17] now states that a basis C of master integrals can be found such 
that the system of differential equations in the kinematic variables Xj is given by 

diC{xj,e) = eAi(xj)C{xj,e), (4.1) 

where di = d/dxi. The C{xj,e) denote the N master integrals and Ai{xj) are N x N 
matrices which are independent of e. It turns out that eq. (4.1) can be expressed in a 
compact form 

dC{xj, e) = e [d A(xj)')C(xj, e), (4.2) 

^The analytic results in section 6 contain only 01 — 04. The results of the “mass-flipped” integrals (see 
section 6 and [ 31 ]) contain also 05 and oe. 
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with the function A determined by the differential diA = Ai. We note that A, together 
with the boundary conditions, completely determines the solution to a master integral. 
The master integrals in such a basis have in turn several pleasant features: First, the 
solution decouples order-by-order in the e-expansion. Second, it is given by pure functions 
to all orders in e. Consequently, assigning a weight —1 to each power of the expansion 
parameter e and multiplying each master integral by an appropriate power of e renders 
the total weight of the master integral to be zero to all orders. Third, the solution can be 
expressed in terms of iterated integrals. If the coefficients Ai{xj) are rational functions of 
the Xj , the Goncharov polylogarithms discussed above represent a suitable class of iterated 
integrals to describe the master integrals. We will refer to such a basis as a canonical basis. 

In the absence of a completely general algorithm for the systematic construction of 
the canonical basis, the procedure of finding such a basis requires a certain amount of 
experience and experimentation. In our case, we start from a “traditional” basis that 
consists of undotted and singly-dotted integrals, and compute them up to terms that involve 
functions of weight two. For this task, alternative approaches like Feynman parameters 
or Mellin-Barnes representations [41, 42] have to be used. Afterwards one plugs these 
expressions into seemingly more complicated integrals like the ones in figures 2 and 3 and 
investigates if the resulting expressions are uniform or even pure. This method is mostly 
based on trial and error, but has proven to be successful as we show below. 

In the case at hand, many master integrals can be adopted from several B —?■ Trvr 
calculations [8-10,43]. In order to describe the yet unknown ones in the canonical basis, a 
set of 39 integrals is needed. We obtain the following expressions for the canonical master 
integrals Ci-^g in terms of the integrals / 1 - 42 , which are defined in figures 2 and 3 {x = 


1 — x). 

Ci{u, z) = uz Ii{u, z), (4.3) 

C 2 {u, z) = u{z — l)z hiu, z), (4.4) 

Csiu, z) = uz hiu, z), (4.5) 

C 4 {u, z) = uz l 4 {u, z), (4.6) 

C^iu, z) = u{z - 1) hiu, z), (4.7) 

Cq{u, z) = u{z - 1) lain, z), (4.8) 

C 7 {z) = e{l - e)zl 7 {z), (4.9) 

Csiu, z) = {u + uz) Isiu, z), (4.10) 

Cg{u, z) = uz(jg{u, z) + 21^{u, z)^ , (4.11) 

Cio{u, z) = €^ (u + uz) Iio{u, z ), (4.12) 

C'ii(u, z) = u{z - 1) z) + 2/io(u, z)^ , (4.13) 

Ci 2 ( 2 :) = e^Ii 2 (^), (4.14) 
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h{u, z) 



Ifi(u, z) 



h{u,z) 




Is{u,z) 



Ii7{u,z) 



hiiii, z) 



Ii4(z) 




h2 



Figure 2: Part I of the basic integrals needed in the construction of the canonical basis: 1,... ,4 
denote the incoming momenta qi,...,q 4 . The double/curly/dashed line represents a propagator 
with mass mf,/mc/0. The dot on a line indicates a squared propagator. 
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1+4 


1 + 4 



hsiu.z) 





h7(u,z) 


1 + 4; 





l2fi[u,z) 


l27{z) 


hs{u, z) 






1 + 2 




1 + 4 


hi(z) 


1 + 4 



h>i(u,z) 


h9(u,z) 


ho(u, z) 



Figure 3: Part II of the basic integrals needed in the construction of the canonical basis. All 
symbols have the same meaning as in figure 2. 


Ci 3 (u, z) = UZ Il 3 (u, z), 

(4.15) 

Cu(z) = e^zluiz), 

(4.16) 

Ci5 = Ii5 , 

(4.17) 

Ciq{u, z) = UZ Ii6(m, z) , 

(4.18) 

Cn{u^ z) = UZ In{u, z), 

(4.19) 
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Cisiu, z) = (1 - uz) (hsiu, z) + -^In{u, z) + z) ) , 

V rnf rnf J 

Ci9(2:) = zlig{z), 

C2oiz) = z (^I2 q{z) + 2/19(2;)^ , 

C21 {u, 2:) = (1 - uz) hi (u, z), 

C22 = I22 , 

C23 {u,z) = uz I23 {u,z), 

C2a{u, z) = uz h^iu, z ), 

C 25 {u, Z) = (1 - uz) [hdiu, z) + -^hiiu, z) + ^hsiu, 2;) ) , 

V w-b / 

C26{u, z) = (1 - uz) heiu, z), 

C 27 {z) = zhriz), 

C28{u,z) = uz hs{u,z), 

029 {u,z) = uz l 29 {u,z), 

C 3 o{u, z) = l-e^ uuz'^ (hi{u, z) + ho{u, z) - -—- ] . hiz)] , 

2 \ e m^uz J 

Csiiz) = zl 32 {z), 

C 32 {z) = ^/z (^l 33 iz) + 2/32(2:)^ , 

C 33 {z) = e^zhhz) , 

(734(2:) = 6^2/35(2), 

C 35 {u,z) = uz he{u,z), 

C3q{u, 2) = (1 - uz)'^ l37{u, 2), 

C'37(u, 2) = (1 - uh hsiu, 2), 

C38 {u,z) = uz I39 (ii, 2), 

(739(«, 2) = |it2 [1 - (1 - uz)p] ho{u, 2) 


(4.20) 

(4.21) 

(4.22) 

(4.23) 

(4.24) 

(4.25) 

(4.26) 

(4.27) 

(4.28) 

(4.29) 

(4.30) 

(4.31) 

(4.32) 

(4.33) 

(4.34) 

(4.35) 

(4.36) 

(4.37) 

(4.38) 

(4.39) 

(4.40) 

(4.41) 


with 

1 — \/{2 — uz)'^ — 42(1 — uO) 
^ 1 — uz 


(4.42) 


Note that the master integrals have to be evaluated to (7(e^) since the two-loop amplitude 
contains poles up to 1/e^ stemming from the infrared and ultraviolet regions. A few 
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exceptions are 6*26,38 and 6*39 which only enter the hard-scattering kernel to order 0 {e^) 
and O(e^), respectively. 


5. Boundary conditions 


Before we present the differential equations, we specify the boundary conditions that are 
used to completely fix the solution. In the simplest cases, the master integrals vanish in 
a specific kinematic point. This is the case for 6*13^38^39, which vanish in n = 0 , whereas 
63,4,5,6,16,17,23,24,28,29,30,35 vanish in u = 1 . Moreover, 619,31,32 vanish in z = 0 , whereas 
67,14,33,34 vanish in z = 1 . In other cases we find special relations between integrals, 
that hold either in general, or in certain kinematic points, and can be used as boundary 
conditions. Examples are the relation 626 = -7”^ 621, or the following relations that hold 
in u = 1 , 

„ U—>-1 „ „ tJ—>1 

t^8 -1 t>l9 , Oio -> Oi9 , 

69 ^ 620 , 611 6^2^ , ( 5 . 1 ) 

where the symbol is used for the corresponding “mass-flipped” integral, in which 

rric o ruft and 53 0 q^, see section 6 for more details. Hence, the integrals 6^20 
easily obtained from 619,20 or from [ 31 ]. Relations that have a similar structure than ( 5 . 1 ) 
hold in z = 1 for 

612 ^ 622 , 627 615 . ( 5 . 2 ) 


For the remaining integrals we either use that they assume simple, closed forms that are 
valid to all orders in the e-expansion, or asymptotic forms as rt —>• 0 or z —>• 0 . Examples 
of the former type are (see below in section 6 for the precise definition of Ci) 


615 

622 


r4(i 


r2(i- 


-e)r(l-4e)r(l + e)r(l + 2e) 
4r(l - 3e)r(l - 2e) 


636 



X 


^ (1 _ -e)T{l + e) 2Fi{l,l + e; 2 -e;u + uz) 

^^:;L\';r^ r(l-e)r(l + e)2Er(M + e;2-e;. + g) 


( 5 . 3 ) 


where for 6*35 we give the result for each loop separately, such that also the boundary 
conditions for 621,37 can be read off. Asymptotic expansions as u —)• 0 or z —0 were 
derived by means of MBasymptotics .m [44] for 


620 _ 1 _ ^ 2C3 6^ - ^ 6^ + 0(6^ .) , 

61 ^ ^ ^ i 

+e^ [^ In^(M) (^ Gi{z) - 60(2;) - i7r)ln(M) ^ iir Go{z) - i 7 rGi{z) -b ^ 60,0(2;) 

o o Z 4 o 
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-I Go, 1 ( 2 ) - 1 G,,„( 2 ) + J Oi,i( 2 ) - 5^1 

+c^ [-j^ln^(u) + (Go( 2 ) - ^Gi(z) + 2i7r)ln^(i£) + (4 27rGi(2) - SiTrGoiz) 
y o 

3 8 377r^ \ 1 / \ 377r^ 

“2 + 2 Go,i(^;) + 2 ^ 1 , 0 ( 2 ) - - Gi^i{z) + ln(M)- Go(z) 

37 yr 2 5 I 

H— t^Gi{z) + -iTrGofi{z) - 3i7rGo,i(2;) - 2 iTrGifi{z) + 4i7rGi^i{z) + - Go,0,0(2) 
lo 4 o 

Go,0,1(2) - ^ Go,i,o(2) + 2 Go,i,i( 2) - ^ Gi,0,0(2) + 2 Gi,o,i(2) + 2 Gi,1,0(2) 

8 ^ / N 17^ 7 , 3 , 

— Gi 1 i(z)-Cs- 

3 ’ ’ ^ ^ 6 ^ 12 ^ 

+e‘^ In^iu) + (y Gi(z) - ^ Go{z) - yvr) 111^(11) + {Qi 7 rGo{z) - 8 i 7 rGi{z) 

+3 Go,0(2) — 4 Go,1(2) — 4 Gi 0(2) + — Gi 1(2;)--— ) In^(tt) + ( —-— Go{z) 

3 9 D 

7471^ 34 

-^ Gi( 2 ;) - 5i7rGo,o(2) + 12i7rGo,i(2) + 8 77rGi,o(2) - 16ivr Gi,i( 2 ;) + y Ca 

“2 ^o,o,o(2) + 6 Go,o,i(2) + 6Go,i,o( 2) — 8Go, 1 , 1 ( 2 :) + 2Gi,o,o(2) — 8Gi, 0 , 1 ( 2 ) 

Qo 7 17 

-8 Gi,1,0(2) + y Gi,1,1(2) 4 - - in^) In(ii) - — Go{z) + 2i7r^ Gi(2) - 8 Gi,0,1,1(2) 

3571^ , , 3771^ n 2 ^ \ 7471^ ^ \ 3 . 

— 7^ Go 0(2) 4 - - — Go,i(2) 4 - 371 Gi 0(2)- — Gi 1(2) 4 - - ^ 7 ^Go 0,0(2) 

24 0 9 4 

-5i7rGo,o,i(2) - 4z7rGo,i,o(2) + 12 ztt Go,i,i(2) - 2 ztt Gi,o,o(2) + 8771 Gi,o,i(2) 

3 13 

+6 77rGl,l,o(2) - 16 771 Gl,1,1(2) + - Go,0,0,o(2) - - Go,o,o,i(2) - - Go,0,1,o(2) 

o 2 2 

1 827r^ 

4 - 6 Go,o,i,i( 2 ) 4 - - Go,i,o,o(2) 4 - 6Go,i,o,i(2) 4 - 6Go,i,i,o(2) — 8Go,i,i,i(2) 4 - 

“2 Gi,o,o,o(2) 4 - 2 Gi,o,o,i(2) 4 - 3 Gi,o,1,0(2) — 8 Gi,i,o,i(2) — 8 Gi,1,1,0(2) 

32 17 34 

+ y Gi, 1 , 1 , 1 ( 2 ) —— Go( 2)C3 4- y Gi(2)C3 — I 2771 C 3 ] 4- 0(e^,u ), (5.4) 


6"3“4 + o[-ll„(,,)-looO)-lo,W-li,l 

+e^ \v?{u) + (^Go(2) + ^Gi(2 ) + y 7 i)ln(n) + -^ 77 iGo( 2 ) + y 7 iGi( 2 ) 

3 6 3 3 12 3 

+e^ ln 3 (n) + (-^ Go(2) “ ^ ^ + (“ J ^ 77 iGi( 2 ) 

Go,o(2) - I Go,i(2) - ? Gi,o(2) - ^ Gi,i( 2) - ^) ln( 77 ) - Go(2) 

Il 7 l 2 112 4 

—y- Gi( 2 ) - y ^71 Go,o( 2 ) - 3 ^^i Go,i( 2 ) + 3 ^^i Gi,o( 2 ) - 3 ^^i Gi,i(2) 
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-^GoM^) - « Go,o,i(^) - I Go,1,0(2) - I Go,i,i(z) + I Gi,o,o( 2) - ^ Gi,o,i(2) 

24 : b 0 6 b 6 

Gi,i,o( 2 ) - - Gi,i,i(z) “ C3 - I 

4 4 18 R 2 8 

+e^ [- In'^(M) + (- Go( 2) + y Gi( 2;) + - m) In^(w) + (- ivr Go(2) + 3 Gi(2) 


+ 3 Go,o( 2 ) + g Go,i( 2 ) + g Gi,o( 2 ) + y Gi,i(z 


yy) ln2(ix) + (yy Go(2) 


22 vr2 4.^,. 8 .^,, 16.^,, 

H—^Gi( 2 :) + -z7rGo,o(2) + -z7rGo,i(2) - -z7rGi,o(2) + yzvr 61 , 1 ( 2 ;) 

1 9 2 R If) R 

+ g Go,0,o(2) + - Go,0,i(2) + - Go,1,o(2) + - Go,l,l(z)-^ Gi,o,o( 2 ) + - Gi,0,i(2) 

+ ^ Gi,i,o( 2 ) + y Gi,1,1(2;) + y Cs + ln(n) + ^ i 7 r^ Go(z) + ^ zvr^ Gi( 2 ;) 

llvr^ „ , , llvr^ „ , . 571^ 22?!^ r \ 4 • ^ \ 

4 —y Go,o(2) + yy Go,i(2)- — Gi^o{z) H- — Gi^i{z) + y Go,0,0(2) 

124 10 8 

+ 3 iTT 60,0,1(2) - -77160,1,0(2) + -77160,1,1(2) + y 77161,0,0(2) - -77161,0,1(2) 

— —771 61 1 0(2) + — 771 61 1 l(z) + — 60,0,0,0(2) + - 60,0,0,1(2) + - 60,0,1,0(2) 

o o 24 o o 

2 5 2 2 8 4 

+ 0 Go,0,1,1(2) — - 60,1,0,0(2) + - 60,1,0,1(2) + - 60,1,1,0(2) + - 60,1,1,1(2) — - 771 Cs 

c5 D o o o o 

19 ^ ^ 10 ^ 7 ^ 8 ^ , , 16 ^ , , 4971^ 

+ y Gi,o,0,o(2) - y 61,0,0,1(2) - - 61,0,1,0(2) + - 61,0,1,1(2) - y 61,1,0,0(2) + yy 

Q Q 09 1 V 94 

+ 0 Gi 10,1(2) + - 6111 0(2) + — 61111(2;) + — 60(2) Cs + “^ Gi( 2;) C3] + G(e®, n) 

3 3 4 6 3 

(5.5) 


6i8 “=V [61(2) 111 ( 77 ) - 60,1(2) + 61,1(2)] 

+e^ [(60,1(2) - 6 6 i,i( 2 ))ln( 77 ) - 61(2) ln^( 77 ) + y 61(2) + 560,1,1(2) - 661,1,1(2)] 

2 271 ^ 

[3 111^(77) 61(2) + (6 61,1(2) — 60,1(2)) ln^(77) + (y“ Gi(2) + 60,0,1(2) 

571 ^ 

-660,1,1(2) - 461,0,1(2) + 2861,1,1(2)) ln(77) + 5C3 Gi( 2) - -y 60,1(2) 

TT^ 

+ y Gi,i(2) — 60,0,0,1(2) + 60,0,1,1(2) + 460,1,0,1(2) — 2260,1,1,1(2) + 461,0,0,1(2) 
~ 4 61,0,1,1 (2) + 2 61,1,0,1(2) + 28 61,1,1,1(2)] + G(e^, 77), ( 5 - 6 ) 

G25 = + (Gi(2) — 60(2)) 111(77) + 60,0(2) — 60,1(2) — 61,0(2) + 61,1(2) + —] 

+e^ [— ln^(77) + (- 60(2) — 3 61(2)) 111^(77) + (3 60,1(2) + 2 61,0(2) — 6 61,1(2) 

-TT^) ln( 77 ) - ^ 60(2) - ^ 61(2) - 360,0,0(2) + 60,1,0(2) + 360,1,1(2) 

2 6 
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+2 Gifl^i{z) + 2 Gi^i^o(^) ~ 6 Gi^i^i(z)] 

+^4 [7 ^ _ 5 ^ (1 -bGo,i{z) - 3Gi,o(^) 

Ttt^ "iTr^ 

+lAGi^i{z) H—^) In^(u) + (Svr^ G*i(z)-^ Gq{z) — Go,o,o(^) + Go,o,i(^) 


— 10 Go,1,1(2;) + 2 Gi,0,0(2;) — 6 Gi,o,i(^) — 4 Gi,1,0(2;) + 28 Gi,i,i(z) + 6 C3) ln(n) 

1 Qtt^ *^7r^ 1 

-6C3 Go{z) + 3 C 3 Gi(z) + — Go,o(^) - — Go,i(^) - — Gi,o(z) + — Gi,i(z) 

5 2 3 3 

+10 Go,o,0,0(2^) — Go,o,o,i(^) — 3 Go,o,i,o(^) + Go,o,1,1(2;) — 2 Go,i,o,o(^) — 2 Go,1,1,0(2;) 

— 10 Go,i,i,i(^) — 2 Gi,o,o,o( 2 ;) + 2 Gi,o,o,i( 2 ;) + 2 Gi,o,i,o( 2 ;) — 6 Gi,0,1,1(2;) 

-I ^ 4 

+2 Gi,1,0,0(2;) — 4 Gi, 1,0,1 (z) — 4 Gi, 1,1,0(2;) + 28 Gi,i,i,i(z) 3 —(^- 7 ) 


6. Results 


In order to facilitate the presentation of the results we write the master integrals as 


G= 


( 6 . 1 ) 


with an integer n that denotes the sum of all propagator powers, such that the integral G is 
dimensionless. Our integration measure is f d^k/{2TT)^ per loop and we use the pre-factor 


Sr 


1 

( 47 r)^/^ r(l - e) ' 


( 6 . 2 ) 


Besides the integrals defined in section 4, the QCD amplitude also contains the same set of 
integrals but with rric •H' mb and -H- 54 . We will refer to these as “mass-flipped” integrals 
and denote them as G'^, see section 5. However, we note here that in order to define G^ 
we factor out an appropriate power of m^, rather than me- 

As stated earlier the QCD amplitude requires terms of order G(e^) for most of the 
integrals. However, in order to keep the paper at a reasonable length, we only give terms 
up to order 0{e^) explicitly below. If desired, terms of weight four can be derived from 
the A and the boundary condition, which we actually give to weight four. Moreover, we 
refrain from presenting the “mass-flipped” integrals explicitly. They can be obtained by 
letting 2 ; —)• 1 /z, keeping in mind that analytic continuation is done via 2 ; —)• 2 ; — fry, with 
infinitesimal ?y > 0. We provide the results to all integrals, including the “mass-flipped” 
ones, to order G(e^) in electonic form in [31]. 

Last but not least, instead of dealing with one large 39 x 39 system of equations, we 
solve each topology separately and therefore deal with several, smaller matrices A, which 
we collect in appendix A. This finally puts us in the position to present the analytic results 
to the Gi_ 39 . 
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6.1 Cl — Ci2 

We start right away with the largest topology, which contains twelve integrals, 

C = {Ci, C 2 , C 3 , C 4 , C 5 , Ce, ^7, Cs, Cg, Cio, Cii, C 12 } . (6.3) 


The corresponding matrix is ^i_i 2 . Taking into account the boundary conditions specified 
in the previous section, the solution to the twelve integrals reads 

Cl + e [-g Go{u) + g Go{z) -^Gi{z) + \ ivr] 

Go{z) Go{u) + ^ Gi{z) Go{u) - i7rGo(u) + ^ iirGoiz) - i 7 rGi{z) 

+ 2 Go{z) Ga2iu) — 2 Cl(^) Ga2{u) + 2 Ga2{u) + - Go,o('w) + g Co,o(^) — 2 Co,l(^) 

- 2 Ci,o(^) + 3 Ci,i(z) - 2 Ga 2 fl{u) - 

+6^ [-3ivr Go{z) Go{u) + divr Gi(z) Go(rt) - ^ Go,o(^) Go(u) + 2 Go,i(z) Go(n) 
+2Gi,o(z) Go(n) - ^ Gi,i(z) Go(u) + ^ Go(ix) - ^ Go(^) + mGo(2) Gi(u) 

+ J^g Gi{z) +2'^^ Go{z) GaaCtt) — 2 iTT Gi( 2;) Gaj (11)-Gaa (u) + 2 Go(2;) Go,0(11) 

Gi(2;) Go,o(u) +4 i7rGo,o(u) + Gi(tt) Go,0(2:) - Ga2(u) Go,o(z) + ^ i 7 rGo,o(z) 

Ga2(u) Co,1(2) - 3 f 7 rGo,i( 2 :) - 2 Go(2;) Go,a2iu) + 2 Gi(z) Go,a2iu) - 2 ivrGo,a2(ii) 

Ca 2 (tt) Gi,o( 2) - 2f7rGl,o(2) + 2Ga2(ll) Gi,i( 2 :) + 4f7rGl,l(z) - Gi( 2 ;) Gl,ai(M) 


+Co(2) Gi,02(11) — Ci(2;) Gi,02(11) + inGi,02(11) — - Go{z) Ga2,o(ii) + 2 Gi(z) Ga2,o(u) 

1 1 18 
-2i7rGo2,o(ii) - 2 Co(2) Go 2,02(11) + 2 C02,02(11) - 2 *^Co2,a2(ii) - 3 Co,0,0(11) 

13 3 1 

+ g Co,o,o(2) — - Co,0,1(2) — - Go,i,o( 2) + 2 Go,1,1(2;) + 2 Go,02,0(11) ~ 2 Ci,0,0(2) 

g 

+2 Gi,0,1(2) + 2 Gi,1,0(2) — 3 Gi,1,1(2;) — Gi,ai,o(ii) ~ Cl,02,0(11) + 2 G02,0,0(11) 


+ 2 Ca2,02,0(11) —0“ ~ 12 ’ 


(6.4) 


C2 ^ ^ *^0(2) “ ^ Ci(2;) - ^ ivr] 


1 


1 


+e [-Go( 2 ) Go(u) + - Gi( 2 ;) Go(ii) + -i 7 rGo(it) + —i 7 rGo( 2 ) + -i 7 rGi( 2 ;) 


^6 


12 


“2 Coi(n) + - Go,0(11) + ^ Co,0(2) + - Co,i(2) + - Gi,o( 2) + - Gi,i(2;) 

1 „ . , llvr^ 

“2 Cai,o(ii) + 

+e^ i7rGo(2) Go(u) - ^ i 7 rGi ( z ) Go(ii) - ^ Go,o(2) Go(ii) - \ Co,i(2) Go(it) 
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-‘^Gifl{z) Go{u) - ^Gi,i(z)Go(u) - - ^^^Gq{z) +i'KGo{z)Gi{u) 

Gi(z) - ^ GaAu) - 3 Go{z) Go,o(«) - 3 G^{z) Go,o(n) - 3 ztt Go,o(n) 

+Gi(u) Go,0(2:) - ^ i7rGo,o(z) - ^ Ga^{u) Go,1(2:) - ^ i7rGo,i(2:) + 2Gi(2) Go,ai(tt) 
G'ai(tt) Gi,o( 2:) + ^ i 7 rGi,o( 2 :) + 2 Gai{u) Gi,i( 2:) - ^ i 7 rGi,i( 2 :) - Gi(z) Gi^aAu) 
+^0(2) Gi,Q,2(n) — Gi( 2:) Gi,a2('u) + Gi,a2('w) — - Go{z) Gai,o{u) + 2 Gi( 2;) Gai,o(^^) 

+ 2 Gliz) Gai,ai{u) — 3 Go,0,o(^i) — ^ Go,0,o(^) — g <^0,0,1(2:) — - Go,1,0(2) 

— 3 Go,1,1(2) + 2 Go,ai,o(^^) + g Gi,o,o(2) — 3 Gi,o,i( 2) — - Gi,1,0(2) — - Gi,1,1(2) 

“Gi,ai,o(lt) — Gi,a2,o('w) + 2 Gai,0,o('w) + - Gai,ai,o(^)- ^ ~ 4 G{e^) , 

( 6 . 5 ) 

2 2 

G3 =e^ [ivrG,2(n) Go(2) - G,2,o(u) Go(2) + 2 Ga„a,{u) Go{z) + y Go(2) + y G,2(n) 

+2 Ga2{u) Go,0(2) — Ga2{u) Go,1(2) — Ga2{u) Gifl{z) — 2 Gi( 2) Ga2,a2(^) 

+2 ivr Ga2,a2iu) + Go,0,o(2) — Go,1,o(2) — 2 Ga2,a2,o('l^) + 2 Cs] + O(e^) ) (6-6) 

G4 =e' [Go(u) Go(2) - Ga2{u) Go(2) - i 7 rGo( 2 ) + Gi( 2 ) Ga2(^x) - mGa2(n) - Go,o( 2 ) 

vr^ 

+Go,i(2) + Ga 2 ,o{u) + —] 

6 

+e^ [4i7rGo(2) Go{u) + 3 Go,0(2) Go(tt) — 4 Go,1(2) Go(n) — 2Gi,o(2) Go('u) 

TT^ 37 r 2 TT^ TT^ 

— y Go(tt) H —— Go{z) — 2 ivrGo(2) Gi(u) + y Gi(tt) —— Gi{z) 

- 3 i 7 rGo{z)Ga 2 {u) + 4 i 7 rGi( 2 ) Ga2(M) + Ga 2 {u) - 4 Go( 2 ) Go,o(n) 

—2 Gi(u) Go,0(2) — 2 Ga2{u) Go,0(2) — Sin Go,0(2) + 2 Gi(tt) Go,1(2) 

+ 3 Ga2(u) Go,i(2) + 4i7rGo,i(2) + 4 Go(2) Go,a2('^) ~ 4 Gi(z) Go,a2('^) + 4 :ivrGo,a2 (u) 
+2 Go(2) Gi,o(u) + 3 Ga2 (u) Gi,o(2) + 2 ivr Gi,o(2) — 4 Ga2 (^i) Gi,i(z) + 2 Gi,o,o(2) 

—2 Go( 2) Gi,a2('w) + 2 Gi( 2;) Gi,a2('u) — 2 ivrGi,a2(w) + 3 Go(2) Ga2,o{u) 

—4 Gi( 2;) Ga2,o('w) + 4 ivr Ga2,o{u) ~ 3 Go(2) Ga2,a2 (u) + 3 Gl(z) Ga2,a2 (u) 

— ?>i'KGa2,a2{u) — 2Go,0,o( 2) + 3Go,0,l(2) + 3 Go,1,0(2) — 4 Go,1,i( 2) — 4Go,a2,o('*^) 

~2 Gi,0,1(2) + 2 Gi,a2,o(^i) ~ 4 Ga2,o,oiu) + 3 Ga2,a2fl{u) ~ Cs + 3 + O(e^), ( 6 . 7 ) 

2 2 

C5 =e^ [Gauoiu) Go( 2) - y Go( 2) + y Ga, {u) + G,, («) Go,l(2) + Ga, (n) Gi,o(2) 

—2 Gi( 2;) Ga^,ai (u) — Go,1,0(2) — 2 Ga^^aifliu) + 2 (^3] + G(e^), (6.8) 
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Ce =e^ l-Go(u) Go(z) + Gi(z) - Go,i(z) + G,,,o(n) - y] 

2 2 

[G’o,o(^) Gq{u) + 4 Gq^i{z) Go{u) + 2 Gi^o(^) G'o(tt) + y Gq{u) + y ^0(2;) 

2 2 2 

—^ Gi{u) + y Gi(z) —— Gai{u) + 4 Go(2:) Go,o(^) ~ 2Gi{u) Go^i{z) 

-GaAu) Go,li^) - 4Gi(z) Go,ai(w) - 2Go(z) Gi,o(n) - Ga^u) Gi,o{z) 

-4GaM Gl,l{^) + 2Gl{z) Gl,aM - Go{z) GaM^) - 4Gl(z) GaM^) 

+3 Gi{z) Gai^ai(u) + Go,0,1(2) + Go,1,0(2) + 4Go,1,1(2) — 4Go,ai,o(^) + 2 Gi,0,1(2) 

+2 Gi,ai,o(^i) ~ 4 Gai,0,o(u) + 3 Gaj,aj,o('u) — ^3] + G(e^), (6-9) 

Gy =e G o( 2) + [—Go,o(2) — 2 Gi,o(2) + y] 

[y Go(z)- — Gi(z) + Go,o,o(2) + 2 Gi,o,o(2) + 4Gi,1,0(2) — 2^3] + G(e^), 

( 6 . 10 ) 


G8=e[-Go(n)-Gi(2)] 

+e^ [4Go(it) Gi( 2 ;) — Gai(u) Gi(z) + 4Go, 0 ( 11 ) + 4Gi,i( 2 ;) — Gai,o('*^) + y] 

+e^ [—16 Gi,i( 2 ;) Go(it)-y Go(it)-y Gi(z) + y Ga^(u) — 16 Gi( 2 ;) Go, 0 ( 11 ) 

+6 Gi(z) Go,ai (u) + 4 Gai (u) Gi,i(z) + 4 Gi( 2 ;) Ga^,o(u) — Gi(z) Ga^,ai (u) — 7(3 
— 16 Go, 0 ,o('*^) + 6 Go,ai,o('*^) ~ 16 Gi, 1 , 1 ( 2 ;) + 4 Ga^,o,o(^) ~ Gai,ai,o(^)] + G(e^) , 

(6.11) 

Gg = — 1 + e [4Go('u) + 4 Gi(2;)] 

+e^ [—16 Go(ii) Gi( 2 ;) + 6 Gq^ (u) Gi(z) — 16Go, 0 ( 11 ) — 16Gi,i(z) + 6Gai,o(?^)-y] 

+e^ [64 Gi,i(2;) Go(u) H- - — Go(it) H- - — Gi( 2 ;) — vr^ Ga^(u) + 64Gi(2;) Go, 0 ( 11 ) 

—24Gi( 2;) Go,ai (u) — 24 Gai(u) Gi,i(2;) — 24Gi(2;) Gai,o{u) + 6 Gi(2;) Gai,ai (u) + 20 (,3 

+64 Go, 0 , 0 ('^^) ~ 24Go,ai,o('^^) + 64 Gi, 1 , 1 ( 2 ;) — 24Gai,o,o('^^) + 6Ga^,aj,o(it)] + G(e^), 

(6.12) 


Gio =e [-Go(n) + Go( 2) - Gi(z) + in] 

+e^ [—2 Go( 2 ) Go(m) + 4 Gi( 2 ;) Go{u) — Ain Go{u) + 2 ivr Go(2) — 4 ivr Gi(2;) 
+Go(2) Ga2{u) — Gi{z) Ga2{u) + inGa2{u) + 4Go,o(ll) — 2Go,i(2) — 2Gi,o(2) 

+4Gi,i(2;) — Ga2,o{u) --—] 
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+e^ [- 8 i 7 rGo{z) Go{u) + l 6 i 7 rGi{z) Goiu) - 4 Go,o(- 2 ) Go{u) + 8 Go^i{z) Go{u) 

TQtt^ fi IQTr^ 

+8 Gifl{z) Gq{u) — 16 Gi^i{z) Gq{u) H-^— Go{u) — -n Go{z) H-^— Gi(z) 

llvr^ 

+2i7rGo{z)Ga2{u) - 4ivr Gi(z) Ga^(u) -^ Gaj(u) + 8 Go(^;) Go,o (u) 

-16Gi(2;) Go,o(^^) + 16 ^ 7 ^Go,o(^^) + 4i7rGo,o(2;) - 2Ga2{u) Go,1(2;) - 8i7rGo,i(2;) 
-6 Go(2 ;) Go,a2('w) + 6 Gi(2 ;) Go,a2('w) - 6ivTGo,a2(w) - 2Ga2{u) Gifi{z) -8iTrGifi{z) 
+ 4 : Ga2(u) Gi^i{z) + 16 ivr Gi,i(2;) — 2 Go(2;) Ga2,o{u) + 4 Gi(2;) Ga2,o{u) + 4 Go,0,0(2;) 
-4:i7rGa2fliu) + Go{z) Ga 2 ,a 2 {u) - Gi{z) Ga 2 ,a 2 {u) + ITT Ga 2 ,a 2 {u) “ 16Go,0,o(^^) 

—4Go,o,i(2;) — 4Go,i,o(2;) + 8 Go, 1 , 1 ( 2 ;) + 6 Go,a 2 ,o('w) — 4Gi,o,o(2;) + 8 Gi, 0 , 1 ( 2 ;) 

+8 Gi, 1 , 0 ( 2 ;) — 16 Gi, 1 , 1 ( 2 ;) + 4 Ga2,o,o{u) — Ga2,a2,oiu) — 7C 3 ~ + G(e^), ( 6 . 13 ) 


Gii — — l + e[4 Gf){u) — 2 Gq{z) + 4 Gi (z) — 4 zvr] 

+e^ [8 Go( 2 ;) Go(n) — 16 Gi( 2 ;) Go(u) + 16 ivr Go{u) — 8 ivr Go( 2 ;) + 16 itt Gi{z) 
-6Go{z) Ga 2 {u) +QGi{z) Ga 2 {u) -6i7rGa2{u) - 16Go,o(w) - 4Go,o(2;) + 8 Go,i( 2 ;) 

IDtt^ 

+ 8 Gi,o( 2 ;) — 16Gi,i(2;) + 6 Ga2,o{u) H- - —] 

+e^ [32i7rGo(2;) Go{u) — 64z7rGi(2;) Go(w) + 16Go,o(2;) Go(tt) — 32Go, 1 ( 2 ;) Go(u) 

—32 Gi^o(z) Go(u) + 64Gi,i(2;) Gq{u) -^— Gq{u) H-^— Go(2;)-^— Gi(z) 

-12 ivr Go ( 2 ;) Ga 2 (u) + 24i7rGi(2;) Ga 2 (u) + llvr^ Ga 2 (u) - 32 Go ( 2 ;) Go,o(u) 
+64Gi(2;) Go,o('w) - 64i7rGo,o('w) - 16 ivr Go, 0 (^ 2 ) + 12Ga2(n) Go,i( 2 ;) + 32 ivr Go, 1 ( 2 ;) 
+24 Go( 2 ;) Go ,02 (u) — 24 Gi( 2 ;) Go ,02 (^) + 24 ivr Go ,02 (u) + 12 Ga 2 (u) Gifi(z) 
+32ivrGi,o(2;) - 24Ga2(u) Gi,i( 2 ;) - 64ivrGi, 1 ( 2 ;) + 12Go(2;) Ga 2 ,o(u) + 64Go,o,o(w) 
-24 Gi(2;) Ga 2 fi{u) + 24ivr Ga 2 ,o(M) - 6 Go( 2 ;) Ga 2 ,a 2 {u) + 6 Gi( 2 ;) Ga 2 ,a 2 {u) 
-6i7rGa2,a2{u) - 8Go,o,o(2) + 16Go,o,i( 2;) + 16Go,i,o(2) - 32Go,i,i(z) 

—24 Go, 02,0 (^^) + I 6 G 1 , 0 , 0 ( 2 ;) — 32 Gi, 0 , 1 ( 2 ;) — 32 Gi, 1 , 0 ( 2 ;) + 64 Gi, 1 , 1 ( 2 ;) 

—24 Ga2,o,oiu) + 6 Ga 2 ,a 2 ,o('*^) + 20^3 + 4ivT^] + G(e^), (6-14) 

2 2 

Gi 2 =1 — eGo( 2 ;) + [Go,o( 2 ;) + —] + Go{z) — Go,o,o( 2 ;)] + G(e^) • (6.15) 


6.2 Gi 3 — Gi 5 

The new integrals in this topology are G 13 - G 15 . However, in order to close the system of 
differential eqnations nine integrals are needed which we order as follows, 

G = {Gi3 , G 5 , Ce, G 7 , Gi4, Gs, Gg, G 15 , G 12 } . (6.16) 
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The corresponding matrix is Ai 3 _i 5 . The solution to the C 13 - Cis reads 


Ci 3 =e3 [Gi{z) Gi,aAu) + - Go, 1 ( 2 ) Gi{u) - ^ Gi(n) - Go{z) Gi,o( 7 x)] + O(e^), 

(6.17) 

Gi 4 [— Go{z) + Go,i,o(^^) ~ 2 (^ 3 ] + G(e^), (6.18) 

Gi5 = -^-e2^-e3 2C3 + G(e^). (6.19) 


6.3 Gi6 — G 22 

This topology has seven integrals, none of which has appeared in previous subsections. 
They are ordered according to 

G = {Gi 6 , Gi 7 , Gi 8 , Gi 9 , G 20 , G 21 , G 22 } . ( 6 . 20 ) 

The corresponding matrix is ^ 16 - 22 - The solution reads 

Gi6 =6^ [Ga, {U) Go,l (Z) - Ga, (u) Go,l{z) + Ga, (u) Gl,l{z) + Ga, {u) Ggl {z) 

~^Gi{z) Gai,o{u) — Gi{z) Gai,ai{u) + Gi(2;) Ga 2 fi{u) — 2Go,0,l(^) + G(i^,i,o(tt) 

— Gax,ai,o{u) + Ga 2 ,l,o(^i) — 2 (( 3 ] + 0{e ‘^), ( 6 - 21 ) 


Gi 7 — [—Gi{z) Gai{u) + Go 1(2;) — G^j o('^^) + 

D 

2 2 2 

[—2 Go 1 (2;) Gi (tt) —;r (^) —T TT Gq i (z) 

6 6b 

+Ga2iu) Go,1(2) + 6 Gai{u) Gi,i(z) — Ga2{u) Gi^i{z) + 2Gi( 2;) Gi,a^(n) — Ga2,l,o('*^) 
+3 Gi(2:) Gaifl{u) — 2 Gi(2:) Gai^ai(u) — Gi{z) Ga2,oiu) + 3 Gq,0,1(2;) — 4Go,1,1(2;) 

“2 Gl,0,1(2;) + 2 Gl,ai,o(^i) + 2 Gai,0,o(^) + Gai,l,o('U') ~ 2 Gai , ai , o { u ) + 3 ^3] + O ( e ^) , 

( 6 . 22 ) 


Gi 8 =e2 [Go(n) Gi(z) - Go,i(z) + Gi,o(n) + Gi,i(z)] 

2 2 

+e^ [—Go, 1 ( 2 ;) Gi(rt) + — Gi(tt) + — Gi(2;) — 2 Gi(z) Go o('it) + Gq{u) Go,i(2 ) 

6 6 

+Gi{z) Go,ai(«) - 4Gi(z) Gi,o(n) - 6Go(n) Gi,i(z) + Gi{z) Gi,a,{u) -2GoMu) 

+5 Go, 1 , 1 ( 2 ;) + Go,ai,o(u) — 2 Gi,o,o(^) — 2 Gi,i,o(u) — 6 Gi, 1 , 1 ( 2 ;) + Gi,q,^,o(?i)] + G(e^), 

(6.23) 


Gi 9 = - e Gl (z) + [4 Ggi (z) - Go,i ( 2 )] 
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^ Gi{z) — Go,o,i(^) + 4Go,i,i(^) + 6Gi^o,i(^) ~ 16 Gi^i^i(z)] + O(e^), (6.24) 

27r^ 

C*2o = — 1 + e4Gi(z) + [6Go,i(2;) — 16Giji(z)- 

[~3~ 4 ” ® ^0,0,i(^) ~ 24 Go,1,1(2) — 24 Gi,0,1(2) + 64 Gi,1,1(2) + 2 (^3] + G(e^), 

(6.25) 


C21 =e[Goiu) + Giiz)] 

4-e^ [—2 Go(u) Gi(z) + Gaj(u) Gi(2) — 2 Go,o(u) — 2 Gi 1(2) + Gai,o('w) — 

0 

271^ 271^ TT^ 

+e^ [4Gi 1(2) Go('u) H—— Go('u) H—— Gi(2) —— Gq^(u) + 4Gi(2) Go o('u) 

4 4 6 

— 2 Gi( 2) Go,ai (^) ~ 2 Ga^ (it) Gi,i(2) — 2 Gi( 2) Gaj,o(lt) + Gi(2) Ga^^ai (u) 

+ 4 Go, 0 ,o(^) ~ 2 Go,ai,o(^) + 4Gi, 1,1(2) — 2 Gai, 0 ,o(^) + G'ai,ai,o('*^) + C3] + C*(f"^) ) 

(6.26) 

G22=l + €^Y^O(e^). (6.27) 

6.4 G 23 — G 27 

This topology has also seven integrals, of which G23 - G27 are new. The entire topology 
reads 

G = |G23, G24, G25, G7, G26, ^27, ^12! • (6.28) 

The corresponding matrix is ^23-27. The solution reads 

2 2 2 

G23 [~Ga2fi{u) Go(z) + — Go(2) + — Gai (u) + — Gaa {u) + Gaj {u) Go,0(2) 

-Ga^iu) Go,i( 2) - GM Gi,o( 2) + G,, (it) Gi,i( 2) + Ga^iu) Gi,i{z) + Gi( 2) GaMu) 
— Gi{z) Gai,ai{u) + Gi{z) Ga2fliu) + Go,0,o(^) ~ Gq, 1,0(2:) + Gaj,o,o(?^) ~ Gai,ai,o(^^) 
+G'a2,o,o('*^) + 2 C3] + G{e^), (6.29) 

C24 =e" [G'o(n) Go( 2) - Gi(2) Ga,(n) + Go,i(2) - G,,,o(n) + 

6 

+e3 [-Go,o( 2) Go(it) - 4 Go,i( 2) Go(n) - 2 Gi,o(2) Go(u) - ^ Go(n) - tt^ Go(2) 

2 2 2 2 

+y Gi(n) - Gi( 2) + y GM - y Ga2{u) - 4 Go( 2 ) Go,o(«) - G,,(n) Go,o(2) 

+2Gi('tl) Go,1(2) + Gai{u) Go,1(2) + Ga2(u) Go,1(2) + 4Gi(2) Go^ai(u) — 2 Ga^,a^,o('u) 
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+ 2 Go{z) Gi^o('*^) + Gax{u) Gifi{z) + Ga2{u) Gifl{z) + 2 >Gai{u) Gi^i{z) — Ga2,0,o('*^) 

— Ga2 (u) Gl^l(z) — 2 Gl{z) Gl,ai (u) + Go{z) Gaifl{u) + 3 Gl{z) Gaifl{u) + 3 Gai,0,o{u) 

— 2 Gi{z) Gai,ai {u) + Go{z) Ga2,o{u) — Gi{z) Ga2fl{u) — Goflfi{z) — Go, 0 ,l(^) 

- 4 Go,i,i(^) + 4 Go,ai,o(^x) - 2 Gi,o,i(^) - 2 Gi,ai,o(n) - Cs] + O(e^), ( 6 . 30 ) 

C25 [Go('u) {Gi{z) — Gq{z)) + Go,o(^*) + Gofl{z) — Gq^i{z) — Gifl{z) + Gi^i{z) + —] 

+e 3 [ 3 Go,i(z) Go{u) + 2 Gi,o(^) Go{u) - 6 Gi,i(z) Go(n) - vr^ Go{u) - y Go(z) 

Gi(u) - ^ Gi(z) + 3 Go(z) Go,o(^/) - 6 Gi{z) Gofl{u) - Gi{u) Go,i{z) 

0 b 

+Gi{z) Go,ai (u) — Go(z) Gi^oiu) + ^1(2;) Gi^ai(u) — 6 Go,o,o(^) ~ 3 G'o,o,o(^) 

+^0,1,0(5^) + 3 Go,1,1(2;) + Go,ai,o('^^) + 2 Gi,0,1(2;) + 2 Gi,1,0(2;) — 6 Gi,1,1(2;) 
+Gi,ai,o('tf)] + G(e^), ( 6 . 31 ) 

G26 =e [Go(u) + Gi(^)] 

+e 2 [-Go(n) Go{z) - 2 Go(u) Gi(z) + Gi(z) G,,(«) - 2 Go,o(«) - Go,i(z) - Gi,o(2) 

vr^ 

“2 Gi,i( 2 ;) + Gai,o(^^) — 

[Go,o(^) Go{u) + 2 Go,1(2;) Go(u) + 2 Gi,o(2;) Go(u) + 4 Gi,i(z) Go(u) H—y Go(u) 

TT^ 27 r^ TT^ 

+“^ Go( 2 ;) H—— Gi( 2 ;) —— Gai (u) + 2 Go( 2 ;) Go,o(u) + 4 Gi( 2 ;) Go q{u) + Gi o,o(^) 
bob 

-Ga, (w) Go,l (z) - 2 Gi (z) Go,a, (u) - Ga, {u) Gi,o(z) - 2 Ga, (u) Gi,i (z) + 2 Gi,o,l (z) 
“Go(z) Ga^,o('u) — 2 Gi(z) Gaj,o('u) + Gi(z) Gai,ai(u) + 4 Go, 0 ,o(^^) + Go, 0 ,l(^) 
+Go,i,o(^) + 2 Go,i,i(z) — 2 Go,ai,o(^) + 2 Gi,i,o(z) + 4 Gi,i,i(z) — 2 Gaj,o,o(?^) 
+Gai,ai,o(^) + C3] + G(e^), ( 6 . 32 ) 

G27 ~ ~ 4 ^ 2 [“Go,o(2) —[— Go(z) + 2 Go,o,o(^) ~ 2 ^3] + G(e^). 

( 6 . 33 ) 


6.5 G 28 — G 32 

This topology has ten integrals, of which the five integrals G 28 “ G 32 are new. They are 
embedded in the topology as follows 

G = |G28, G29, G30, G7, G26, G21, G22, G31, G32, G12I . ( 6 . 34 ) 

The corresponding matrix is A 28 - 32 . The solution reads 

2 2 2 

G28 = 6 ^ [-Gi,o(w) Go(z) + ^ Go(z) - ^ Gi(n) - ^ Gai(n) - Gi(n) Go,i(z) 

b b b 
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— Gi{z) Go,ai(^) + ^1(2;) Gi^ai{u) + Gi{z) Gai,ai{u) — G'o,ai,o('*^) + 
+Gai,ai,o(^‘)] + O(e^) ) 


( 6 . 35 ) 


2 2 

C 29 =e^ [Go,i {z) Gi (u) + ^Gi(u)-^ Ga, (u) - Ga, (u) Go,i(z) + Gi (z) Go,a, (u) 

0 6 

+Go(z) Gi,o(u) - Ga,(n} Gi,o(z) - Gi(z) Gi,aAu} - Go(z) Ga,,o(u) + Gi(z) Ga,,a^(u) 
~l~Go,l,o(z) + Go,ai,o(u) — Gi,a^,o(u) + Ga^,ai,o(u) — 2 ^ 3 ] + O(e^), (6.36) 

C 30 =e^ l-Go(u) Go(z) + Gi(z) Ga,(u) - Go,i(z) + Ga,,o(u) - y] 

+e^ [2 Go,o('^^) <-^ 0 ( 2 ;) — 3 Gi^o(^) ^ 0 ( 2 ;) — 2 Ga^^o('*^) ^ 0 ( 2 ;) + 2 Ga 3 ,o(u) Go(z} 

jr2 7^2 277-2 27r2 

+2Ga4,o('^^) <^ 0 ( 2 ;) + — Go(2;)- — Gl{u) + — Gl{z) - — Ga^{u) H-^ Gasiu) 

+ “^ Gaiiu) — 2 Gas ('*^) G-l,Q{y/z) + 2 Ga^iu) G-i,o{y/z) + Go{u) Go,q{z) — Go,ai,o(u} 
+2 Go(u) Go, 1 ( 2 ;) — 3 Gi(tt) Go, 1 ( 2 ;) — 2 G^^ (u) Go, 1 ( 2 ;) + 2 G^g (tt) Go, 1 ( 2 ;) + Go, 1 , 0 ( 2 ) 
+ 2 Ga 4 (u) Go, 1 ( 2 ) — Gi( 2 ;) Go,ai(^i) + 2Gas(u) Gi,o(\/z) — 2Ga4(u) Gi,o(\/z) 

+2 Go(u) Gi,o(z) — 2 Gas (^) ^* 1 , 0 ( 2 ) + Gag (u) Gi,o(z) + Ga^iu) Gi,o(z) + Go, 0 , 1 ( 2 ) 

— 2 Gas (u) Gi,i(z) + 3 Gi( 2 ;) Gi,aj (u) — 2 Gi(z) Gas,o(u) + 4 Gi( 2 ;) Gai,ai (u) + Cs 

— 4 Gi(2;) Gas,ai (u) — 4 Gi( 2;) Ga 4 ,ai (u) + 2 Go, 1 , 1 ( 2 ) + 2 Gl, 0 , 1 ( 2 ) — 4 Ga 4 ,ai,oiu) 

+3 Gi,ai,o(^^) ~ 2 Gas, 0 ,o(u) + 4 Gaj,ag,o('u) — 4 Ga 3 ,ai,o(^)] + ^(c^) 7 (6.37) 

G 31 ~ ~ 2 ■*" [~G_i,_i,o(\/ 2 ) + G-i,i,o(\/z) — - Go, 1 , 0 ( 2 ) + Gi,-i,o(\/z) 

+ 2 Gl, 0 , 0 ( 2 ) — Gi,i,o(\/z) + 2 Gl, 1 , 0 ( 2 )] + O(e^), (6.38) 

G 32 =e^ [G_i,o(-v/2) — Gi,o(\/ 2 )] 

+e^ [—2 G_i,_i,o(-v/ 2 ) — 2 G_i,o,o(-\/2) — 4 G-i,i,o(v^) — 2 Go,-i,o(-y/2) 

+2 Go,i,o(vG) + 4:Gi,_i,o(\/ 2) + 2 Gi,o,o(\/2) + 2 Gi,i,o(\/2)] + G(e^) ■ (6.39) 


6.6 G33 and G34 

This topology has seven integrals, of which only G33 - G34 have not yet appeared in the 
previous topologies. The integrals are ordered as 

G = {G33,G34,G7,G22,G31,G32,G12} . ( 6 . 40 ) 

The corresponding matrix is ^33,34. The solution to O(e^) is very short 

G33 =e 3 [Go,i,o(2) - 2 C3] + G(e"), ( 6 . 41 ) 
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(6.42) 


C 34 —— Go(z) + O(e^) • 
b 

At order O(e^) the solution requires also Goncharov polylogarithms of argument -^/z. 

6.7 C 35 

This topology has three integrals, and only C 35 is new. The integrals are ordered as 

C = {C35,Ci9,C2o} . (6.43) 

The corresponding matrix is A 35 . The solution reads 

(535 =e ^ Go(u) + [-2 Go(u} Gi(z) + Ga^) Gi(z} - ^ Go,o(n) - Go,i(z) - ^ Gi,o(u) 

TT^ 

+Gaifl{u) - Yj] 

+e^ [ 8 G 1 i(z) Go{u) + —— G'o(n) + — Gi{u) + — Gi(z) —— Ga^iu) + - Goo o{u) 

12 12 3 6 2 

+ 6 Gi{z) Go,o(u) + Gi(u) Go,i(z) — Ga^ (u) Go, 1 ( 2 ) — 3 Gi( 2 ;) Go,ai (u) — 3 Go,ai,o(^) 
+2 Gi (z) Gi,o(u) - 4 Ga, (u) Gi,i (z) - Gi (z) Gi,„, (u) - 4 Gi (z) Ga, , 0 (^^) + 4 Gi,o,i (z) 

+2 Gi(z) Gai^ai (u) — Go,o,i(^) + 2 (^0,1,o('*^) + 4Go,i,i(z) + - Gi,o,o('*^) + ^ (^ 1 , 1 ,o(^^) 
“G'l,ai,o('^^) ~ 3 Ga^,o,o('^^) ~ Gai,l,o('w) + 2 Gai,ai, 0 ('*^) + Cs] + G(€^) . (6.44) 


6.8 G36 and G37 

This topology has four integrals, of which G36 and G37 are new. The integrals are 

(5 = {(536,(537,(526,Gi2} . (6.45) 

The corresponding matrix is ^30^37. The solution reads 

G 36 =€' [-Go(^) Go(n) + Gi(z) Go(u) - i7rGo(u) + Gi(u) Gi(z) - i7rGi(z) 

~l~Go,i(u) — Go,i(z) + Gi,o(u) — Gi^o(z^) + 2Gi,i(2;)] 

+e^ [divr Gi(2;) Go{u) + Go,o(2;) Go(u) + 2 Go, 1 ( 2 ) Go(u) + 3 Gifl{z) Gq{u) 

“6 Gi i(^) Gq{u) H—— Gq{u) + — Go{z) —— Gi(m) + 2 in Gi(tt) Gi(2:) + — Gi(2;) 

3 6 6 2 

—2 zvr Gi(z) Gai (u) + 2 Go (2) Go,o(^) ~ 2 Gi( 2 ;) Go,o(u) + 2 in Go,o('^) + Go(z) Go,i(u) 
—4 Gi(2;) Go,i(u) + 2ivr Go,i(u) + Gi(u) Go, 1(2) — 2 G^j (u) Go, 1(2) — Go{z) Go,ai (u) 
+Gi{z) Go,ai(u) - inGo^aiiu) + Go(2:) Gi,o(tt) - 4Gi(z) Gi,o(n) + 2z7rGi,o(n) 

+Gi(u) Gi,o( 2 ;) - 2Gai(u) Gi,o( 2 ;) -2Gi{z) Gi,i(m) - 6Gi(n) Gi,i( 2 ;) + 6z7rGi,i(2;) 
+4 Gaj (m) Gi,i(2:) + Gi (z) Gi,ai (u) — 2 Go( 2 ;) Gai,o(^^) + 2 Gi(2;) Gai,o{u) + 3 Gi,o,i(^;) 
—2 in Ga^,o(u) + 2 Gi( 2 ;) Gaj,i(u) — 2 Go,o,i(^^) + Go,o,i(^) — 2 Go,i,o(^^) — 2 Gi,o,i(tt) 


- 22 - 


—2 + 2 Go, 1 , 1 ( 2 ;) + Go,ai,i(?^) — 2 Gi,o,o('it) + Gi,o,o(2:) + Go,i,o(2;) — 2 Gi,i,o(m) 

+ 4 Gi,1,0(2;) — 12 Gi,i,i(z) + Gi,ai,o(^^) + 2 Gq,^,o,i(u) + 2 Gq,^,i,o(u) + - zvr^] + G(e^), 

( 6 . 46 ) 

G37 =e [-Go(^) + Gi{u) + Gi(z) - ITT] 

[Go(2) Gi(11) — 2 Gi(z) Gi(it) + 2ZTT Gi{u) + 2zvr Gi(z) — Go(2;) Gai(u) + Gai,i(w) 

271 ^ 

+Gi(2;) Gai (u) — in Ga^ (u) + Go,0(2;) + Gi,o(2;) — 2 Gi,i(tt) — 2 Gi,i(2;) H — 

+£=" [-2 Gl,l(«) Go{z) + 2 Gl,ai(«) Go{z) + Ga„l(n) Go{z) - Ga,,aM Go{z) 

TT^ 

-—Go{z) - n^ Gi{u) - 4i7rGi(n) Gi(z) - Gi(z) + 2i7r Gi(z) Ga^u) + 4 Gi,i,i(z) 

O 

+ —^ Gai (u) — Gi{u) Go,0(2;) + Gai (u) Go,0(2;) — 2 Gi(u) Gi,o(2;) + Gai {u) Gifi{z) 

+ 4 :Gi{z) Gi,i(u) - 4i7rGi,i(ii) + 4 Gi(it) Gi,i(2;) - 2 Gaj^{u) Gi^i{z) - AinGi^i{z) 

—2 Gi( 2;) Gi,aj (tt) + 2 in Gi,ai (u) — 2 Gi{z) Gaj,i(tt) + 2 i7r Gai,i(u) + Gi(2;) Ga^^ai (u) 
—in Gai,ai{u) — Go,0,0(^) — Gi,o,o(2 ;) — 2 Gi,i,o(2;) + 4 Gi,1,1(11) — 2 Gi,aj,i(it) 

“2 Gaj,i,i(ii) + Gai,ai,i{u) + 2 ((3] + G(e^). ( 6 - 47 ) 


6.9 G 38 and G 39 

These integrals arise from diagrams with a massive quark loop inside a gluon propagator. 
They appeared in a slightly different version already in the calculation of the two-loop tree 
amplitudes in i? —>■ vrvr [9,10], and analytic results were recently derived in [11] as M 28 , 29 - H 
turns out that the results of Gss^sg can be obtained from the latter reference if one adjusts 
the kinematics to the present problem. To be precise, one has to replace 


u —>■ u (1 — z) 


(6.48) 


in the expressions in [11]. That is, in the definition of the canonical basis (cf. (3.30) and 
(3.31) of [11] and (4.40), (4.41) of the present article), and also in the solution, eqs. (4.64) 
and (4.65) of [11]. In particular, the kinematic variable p changes to (z = 1 — 2 ;) 


1 — \/(2 — uzY — 4z(l — uz) 
1 — uz 


(6.49) 


7. Checks 

In order to validate the analytic results presented above, we performed several checks of 
analytic and numeric nature. Those integrals that possess a closed form in terms of hyper¬ 
geometric functions were analytically expanded in e using HypExp [45,46]. Subsequently, 
we re-wrote the resulting polylogarithms and HPLs in terms of Goncharov polylogarithms 
and compared to the results obtained by the differential equation method. 
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For the numerical checks we used a dozen points in the u — z plane. We first evaluated 
the Goncharov polylogarithms that appear in our analytic results numerically with the 
GiNaC-library [47,48]. We also derived Mellin-Barnes (MB) representations, partially 
using the AMBRE-package [49]. The analytic continuation to e = 0 and subsequent numerical 
integration was carried out by MB.m [50]. This worked for almost all cases, even in the 
presence of kinematic thresholds, and yielded agreement to the GiNaC results to 5T0“^^ or 
better. There are, however, a few cases in which the Monte-Garlo integration implemented 
in MB.m failed due to highly oscillating integrands, notably for the integrals C' 28 - 30 ) and 
their “mass-flipped” counterparts (where rric -H- and qs Qi)- In these cases, we 
relied on the sector decomposition method implemented in SecDec [51,52], which yielded 
agreement with GiNaC at the level of 8-10“^ for the highest e-coefRcients in C' 28 - 30 ) and at 
the level of 6 • 10““^ for the highest e-coefficients of their “mass-flipped” counterparts. The 
agreement is several orders of magnitude better for the lower coefficients in the e-expansion. 

Another important point to mention is the fact that the GiNaC results were obtained 
in the canonical basis, whereas most of the MB representations and the SecDec results 
were derived in an “ordinary” basis of un-dotted and singly-dotted master integrals. The 
change of basis was then performed using the Laporta reduction. Having calculated the 
numerics in two different integral bases constitutes another non-trivial check of our results. 

8 . Conclusion 

We obtained analytic results to all two-loop master integrals that are necessary for the 
description of the non-leptonic decay B —)> Dtt at NNLO in QCD factorisation. They are 
expressed in terms of Goncharov polylogarithms of argument u and weights that are either 
integer numbers (0 or ±1), or contain the second kinematic varible, c. It is remarkable that 
six c-dependent weights are sufficient for writing down the entire set of solutions, including 
the “mass-flipped” integrals. 

With the master integrals at hand, the bare two-loop part of the hard-scattering kernels 
Tij{u) in (1.1) is complete. The remaining task consists of renormalising the ultraviolet 
divergences and subtracting infrared divergences via matching from QGD onto soft-collinear 
effective theory. Steps towards this goal are outlined in [13]. Having the hard-scattering 
kernels Tij{u) written in terms of iterated integrals is an optimal choice for carrying out 
the convolution integral with the pion LGDA in (1.1), and it might be feasible to obtain 
the NNLO topological tree amplitude in analytic form. In any case our results constitute 
an important step towards the phenomenology of B —)• Dtt decays at NNLO in QGD 
factorisation. 

Let us compare the integrals in the present work to those recently obtained in [11] 
during the evaluation of the two-loop penguin amplitude. Both are two-loop problems 
with scales u and z. The present integrals are a bit less involved compared to those in 
[11], in a sense that the linear combinations that form a canonical master integral are 
shorter, the occurring weights are fewer, and the choice of kinematic invariants is less 
complicated. The main reason for this is that in the present work the external kinematics 
of the final state contains also the second internal mass, notably rric- On the other hand. 
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the only five-line integral in [11], a two-point function {M22), is in fact a one-scale integral, 
whereas here we encountered several five-line integrals with four external legs which are 
genuine two-scale functions. Moreover, most of our integrals are needed to order O(e^), 
whereas in [11] all but two integrals were required only to order O(e^). 

On more general grounds, it will be interesting to investigate how the canonical basis 
depends on the number of loops, legs, scales, space-time dimensions, and on the external 
kinematics. Every example therefore sharpens our understanding of the patterns that 
such bases follow, with the goal of eventually developing an algorithm for their automated 
construction. 
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A. The matrices A 

Here we list the matrices A for the different topologies. Their entries can all be expressed 
in terms of the following nine logarithms, 

Li =ln(u) , 

L2 =ln(l - u) , 

L3 =ln( 2 :) , 

L4 = ln(l — z) , 

L5 =ln(l - u{l - z)) , 

The matrices A now assume the following compact form, 

- 4 ii - L4 3L3 - 3L4 -2L2 -f-f+Le -L2 + ^-L 4 +f L2 + f-f 

—3L4 —ALi — L3 — L4 L2 — —L2 — L4 —2L2 ~ ^ + is 

0 0 2L2 — L3 ~\~ 2L4 — 2Lq L3 — Lq 0 

0 0 — 2Li — 21^4 + 2L3 —ALi + 2L2 — 2L4 + L3 0 

0 0 0 0 2L2 — L3 2L4 — 2L3 

0 0 0 0 — 2Li — 2X4 + 2L3 

0 0 0 0 0 

0 0 0 0 0 

0 0 0 0 0 

0 0 0 0 0 

0 0 0 0 0 

0 0 0 0 0 



Lq = ln{z + U[l - z)) , 

L7 = In (1 — U (1 — a/z)) , 
Ls =ln (1 - U (1 -F ^/z)) , 
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— is + is — i4 

L3 Lq 

4 4 

0 

L3 L4 

4 4 

3 Lq 3-^3 

2 2 

-^3 4_ -^"4 1 -^6 

2 T 4 T S 

Lq L3 \ 

4 4 

— is — is 4" ^ 

U 

4 

Sis 

2 

L4 1 L5 

4 f 4 

0 

La 

4 

Lb 

4 

0 

2 

0 

0 

— Sig 

La 

2 

Lq 

2 

0 

Li — ^ 

1^1 2 

0 

0 

is 

La 

2 

Lq 

2 

-is 

ia ia 

2 2 

2L3— 2is 

L3 

2 2 

0 

0 

2 2 

4“ 2is — is — 2i4 + is 

r, La Ls 

^ 4- -5- 

is - is 

1/5 1/3 

2 2 

0 

0 

1/5 1/3 

2 2 

0 

-2i4 

0 

0 

0 

0 

4^3 

0 

0 

is 

Li + 4^4 

0 

0 

0 

0 

0 

-eis 

1 

1 

0 

0 

0 

0 

0 

0 

0 

ie — Sis 

il — is 4- L4 

0 

0 

0 

0 

0 

eis - 6ig 

—44 ji + 24^3 — 44^4 

0 

0 

0 

0 

0 

0 

0 

-is / 


(A.2) 


/—2^4 L^ — Li 

0 2X2 — + 2X4 — 2L5 

0 


^13-15 — 


L2 

-L5 


L3 

2 


L5 

2 


—2Li — 2L4 + 2L^ —4Li + 2L2 — — 2L4 + ~Li ~ ^ ^ 


L 3 — L 4 
0 
0 


- 2 L 4 

_ 

2 

0 

0 

0 

0 


0 

2L4 — Ls 
0 
0 
0 
0 


2 I 16-22 — 


/ 2L2 — 2X3 + 2L4 + is is “ is 

— 2i2 + 2i3 — 2i4 + is —2i2 + is — 2i4 + is 

— 2il — 3i2 + 2i3 — 5i4 —il + is — i4 

0 0 

0 0 

0 0 

\ 0 0 


2i4 — 2i3 
Sis “ Sis 
is - is 
0 
0 

is 

-6is 

0 

0 


1/4 _ 

2 2 

^3 _ 

2 2 

^5 ^3 

2 2 

0 

0 

— 4 ivi — 4^4 
0 
0 


2L3 -2L4 0 \ 

ia _ i 

2 2 


—is + is + ig 
L 3 + L 3 — ig 
—2ii — 2i2 + is — 4i4 
0 
0 
0 
0 


0 

0 

0 

—Sis 

0 

0 

0 

0 


L5 L3 
2 2 
is 

_ii 
2 

0 
0 
0 

-is / 

(A.3) 


2i5 — 2L3 

is + 2 i 5 


is 

2 

ia 

2 


— 3il — 3L2 3L3 — 6i4 —^ — is + ^ — 


ia 

2 

ia 

2 

ia 

2 


Sia 

2 


is 

-6is 

0 

0 


is 

-4i4 

0 

0 


0 

-is 
i2 4“ is 
0 
0 

— 2 L /1 — 2 _Zv4 “h -Zv5 

0 


^5 

2 

2 


^3 

2 

L5 

2 


_ _ T J- _ 31/4 

2 -^2+2 2 

0 

0 

4^1 + 4^4 

0 


(A.4) 
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‘IL/2 — 3 X 3 + 2^4 + L<^ ”^3 + + Lq 

— 2 Li — L3 — 2L4 + L5 — 4 Li + 2L2 — 2L3 — 2L4 + L5 —is + — Lq 


5X4 

1 

1 

— 4X1 4 “ Lq — 4X4 



0 

0 



0 

0 



0 

0 



0 

0 


2 2 

0 

2X5 — 2 Lq 

L5 Lz 

2 2 

1 Ls L5 
^2 2 

Lq — Lq 

2X5 — 2X3 

^5 ^3 

2 2 

i '3 J_ •^'4 

2 2 

Xl — X3 + X4 

—6X1 + 6X3 — 6X4 

3 -Li 1 3L3 

2 2 

-2L4 

0 

0 

Lq 

0 

— 2X1 — Lq — 2X4 + Lq 
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